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and about symmetry breaking inflicted on a bilayer by the underlying substrate. 
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I. INTRODUCTION 



A similarity between electrons in monolayer graphene 
and relativistic massless particles has been broadly dis- 
cussed in the literature [311S00. The combi- 
nation of a sublattice composition of electronic Bloch 
states (treated as an "isospin" ) in a single atomic sheet 
of graphite with a linear dispersion in the vicinity of cor- 
ners of the Brillouin zone makes them chiral, like Dirac 
fermions. Experimentally, the chiral nature of charge 
carriers has been deduced from a peculiar sequencing of 
plateaus in the quantum Hall effect [3 |9j [10], while 
the linear dispersion relation has been observed directly 
by angle-resolved photoelectron spectroscopy (ARPES) 

[in[I2|]]am|T3[ig|T3[IB]- ARPES has al- 
ready been used to provide information about the form 
of the dispersion curves, renormalization of spectra by 
electron-electron and electron-phonon interactions, and 
information about quasiparticle lifetimes in the mate- 
rial PH QU EQl EEJ [22] . Here, we point out that 
constant-energy angular maps of photoemission reflect 
the chirality of electrons in graphene. For monolay- 
ers, we show that the recently published ARPES data 
provides evidence for the chirality of carriers in this 
material. We demonstrate that the anisotropy of the 
constant-energy maps may be used to extract informa- 
tion about the magnitude and sign of interlayer coupling 
parameters in bilayer graphene and about the types of 
symmetry-breaking effects produced by the underlying 
substrate or doping. In particular, we demonstrate that 
one can distinguish between two effects that may gener- 
ate a gap in the bilayer spectrum: interlayer asymmetry 
[T31 123 1211 [231 [251 [251 [23 [301 [311 [32] and symmetry 
breaking in the bottom layer resting on a SiC substrate 

H3HSH33J. 
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Whereas the chirality of a relativistic particle is defined 
by its spin, chirality in graphene refers to the sublattice 
composition of plane- wave states of Bloch electrons. The 
honeycomb lattice of monolayer graphene, Fig. 1(a), has 
two sites in the unit cell, labeled as A and B, and there 
is a degeneracy point at each of two inequivalent corners 
K± = ±(47r/3a, 0) of the hexagonal Brillouin zone, also 
referred to as valleys, Fig. 1(b). Near the center of valley 
K + in monolayer graphene, electrons are described by a 
Dirac-like Hamiltonian, 

Hi w hva.q, (1) 

which determines the linear dispersion e = hvq of elec- 
trons in the conduction band and e = —hvq in the va- 
lence band pQ. The eigenstates 'J, within a single valley, 
have different amplitudes on adjacent A and B sites, and, 
following the example of relativistic physics, they may 
be written as a two-component "spinor" = (V'AjV'b)- 
The chirality of a relativistic particle is right-handed if 
its spin points in the same direction as its momentum, 
left-handed if its spin points in the opposite direction. 
By analogy, the relative phase (f> between the wave func- 
tion on sublattice atoms indicates the 'isospin' vector 
a = (cos 0, sin 0, 0) of the chiral state * = (e"^/ 2 , e^/ 2 ) 
of quasiparticles in graphene. ARPES [34j [35] is exactly 
the tool to visualize this state through the angular de- 
pendence of the emitted photoelectron flux. 

The proposed analysis is based upon the standard 
theory of angle- resolved photoemission |34| 135] . In an 
ARPES experiment, incident photons with energy huj 
produce photoelectrons whose intensity / is measured 
in a known direction as a function of kinetic energy 
E P « ft 2 (|P||| 2 +p 2 z )/2m |3J| E5J: Hco = E p + A - e q 
where A is the work function and e q is the energy of 
Bloch electrons in graphene. Conservation of momen- 
tum ensures that the component of the momentum par- 
allel to the graphene surface hp\\ = h(p xi p y ) is equal to 
the quasi-momcntum KK.± + fiq of Bloch electrons near 
valley K±: 

P||=K±+q+G, (2) 
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where a reciprocal lattice vector G = mibx + TO2b2 is 
written in terms of primitive reciprocal lattice vectors 
h>! = (2n/a,2n/y/3a), b 2 = (2n/a, — 2w/V3a), and inte- 
gers mi, 777*2 ■ 

As graphcnc has two incquivalcnt atomic sites, the an- 
gular dependence of the intensity may be accounted for 
by considering two-source interference (a la Young's dou- 
ble slits). Outside the sample, at a position Ro rela- 
tive to the midpoint of the two sources, electronic waves 
e [i P .(R +u/2)-^/2] and e [i P .(R o -u/2)+i0/2] from adjacent 

A and B sites combine. This yields the intensity / of the 
two-source interference pattern, 



cos 



p || .u 



(3) 



where u = (0, a/v3) is the separation of the adja- 
cent sites, and near each corner of the Brillouin zone 
P||.u w 2-k(vti\ — ?Ti2)/3. The first term in the argument 
of Eq. ^ is a phase difference due to the different path 
lengths of electron waves emitted from two sublattices, 
while the second term, —cfi/2, arises from the relative 
phase of the electronic Bloch states on A and B sublat- 
tices determined by the quasiparticle chirality. 
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Table 1: Manifestation of electronic chirality in the 
anisotropy of ARPES constant-energy maps in monolayer 
graphene in the valley K + , where angle ip specifies the 
direction measured from the center of the valley. 

Electrons in the conduction and valence bands at the 
valley K + , determined by the Dirac Hamiltonian, Eq. (JT|), 
differ by the projection of their isospin onto the direc- 
tion of their momentum q = (g cos tp, q sin ip) as described 
by the chiral operator a.xi\ where ni(q) = (cos ip, sin (p): 
(T.ri! = 1 in the conduction band, a.xii = — 1 in the va- 
lence band, as listed in Table 1. Note that, the first term 
in the argument of Eq. pL arising from the path dif- 
ference of electron waves emitted from two sublattices, 
accounts for the relative rotation in the interference pat- 
tern around the six corners of the hexagonal Brillouin 
zone. Figure 1(c) shows a typical calculated dependence 
of the intensity of the photocmission from states (here, 
at energy 1.45eV below the charge-neutrality point [31)] ) 
plotted as a function of wave vector py, in agreement 
with the qualitative prediction of the two-source inter- 
ference picture Eqs. 
numerical data, 



( 1|3 ), summarized in Table 1. The 
Fig. 1(c), appear to be consistent with 
experimentally measured constant-energy maps 

So far, we have discussed the angular dependence of 
the interference patterns, neglecting the effect of trig- 
onal warping. It leads to a triangular deformation of 
iso-energetic lines in the band structure of graphene and 
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FIG. 1: (a) schematic of the monolayer lattice containing two 
sites in the unit cell: A (white circles) and B (black), (b) 
Schematic of the hexagonal Brillouin zone with two inequiv- 
alent valleys K± and the low energy bands e w ±hvq near 
the K + point obtained by taking into account intralayer hop- 
ping with velocity v. Shading indicates the region of occupied 
states up to the Fermi energy 6f , dashed line indicates a typ- 
ical energy of states contributing to photoemission, whereby 
incoming photons of energy hui produce photoelectrons of ki- 
netic energy E p . (c) The intensity of photoemission from 
states at a constant energy 1.45eV below the charge-neutrality 
point |36| in monolayer graphene, plotted as a function of pho- 
toelectron wave vector = (p x ,p y ) parallel to the surface 
of graphene for py covering several Brillouin zones (top) and 
plotted as a function of photoelectron wave vector q = (q x ,q y ) 
in the vicinity of valley K+ (bottom) [note that the origin and 
scale of py and q are different]. Here we use parameter val- 
ues 70 = 3.0eV, so = 0.129, A = 0, and the energy width 
T = 0.24eV. 



e g (— q) 7^ £q(q) asymmetry of the electron dispersion 
around each valley which becomes more pronounced for 
states further from the charge-neutrality point. Another 
perturbation of chiral particles in graphene may be asym- 
metry A = ea — eg of on-site lattice energies , e# due 
to the presence of a substrate, leading to a gap A in the 
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spectrum at low energies. The presence of such asym- 
metry in graphene grown epitaxially on SiC substrate, 
and the possibility of observing its effect within spec- 
troscopic accuracy, has recently been discussed following 
experimental ARPES measurements of the low-energy 
band structure [15j [16]. In Section |TT| below, we show 
that opening an AB asymmetry gap in the monolayer 
spectrum is accompanied by the loss of the chirality- 
related anisotropy of ARPES angular maps at low en- 
ergies, which can be used as an additional test for the 
symmetry-breaking effect induced by a SiC substrate. 

In Section [Tn] we offer a detailed analysis of the angle- 
dependent maps of ARPES of bilayer graphene, which 
is the main goal of this work. First, we analyze angular 
photoemission maps of an ideal "pristine" bilayer, taking 
into account intricate details of its band structure using a 
tight-binding model that employs the Slonczewski- Weiss- 
McClure parametrization of relevant couplings [2j [3] . In 
Sections |III A| and |III B| wc show that angular maps can 
be used to determine not only the magnitude, but also 
signs of the interlayer coupling constants used in the 
tight-binding model. If measured experimentally, the lat- 
ter information may also prove to be useful for general 



studies of bulk graphite. In Section IIIC we analyze the 
influence of inter- and intra-layer symmetry breaking in 
bilayers, and we show that the effect of the interlayer 
charge transfer upon doping can be, in principle, dis- 
tinguished from crystalline asymmetry induced by a SiC 
substrate. 



II. PHOTOEMISSION FROM MONOLAYER 
GRAPHENE 

To produce a quantitative prediction of the photoe- 
mission intensity, we use Fermi's Golden Rule to cal- 
culate the probability of a photo-stimulated transition 
from an initial band state with 2D quasi-momentum 
fik = hK.± + hq and energy e q in graphene to a con- 
tinuum state with momentum hp and energy E p in vac- 
uum [34] . The initial state wave function in graphene is 
written as a linear combination of Bloch wave functions 
on the A and B sublattices with coefficients ipA and ipB, 
respectively: 



^^ e ^R 3 $ (r _ R .) 



* k (r)= ]T ^(k) 

0=A,B 

where R^, R-s are the positions of A and B type atoms, 
and ^(r) is a p z atomic orbital. Then, the intensity / 
of photoemission from states in a given band may be 
written as 



-iG.Tj 



6(E p + A-e q -to), (4) 



alleltothe surface is conserved, q = py— K±— G, Eq. pi. 
The summation with respect to index j = {A, B} takes 
into account coefficients ipA and ips located at atomic po- 
sitions defined by basis vectors ta and tb within a given 
unit cell. The Dirac delta function, containing the work 
function of graphene A, expresses energy conservation. 
Note that, in this paper, we do not model dynamical ef- 
fects that lead to energy broadening jTTJ [T51 [2U1 122| , 
but introduce a Lorentzian 8(. . .) ~ r/(7r[(. . . ) 2 + r 2 ]) in 
the figures with parameter T representing finite energy 
broadening. 

A standard form [371 1SH] of a tight-binding monolayer 
Hamiltonian Hi and overlap-integral matrix Si (that 
takes into account non-orthogonality of orbitals on ad- 
jacent atomic sites), Hi^f = CqSi^, is 



Hi 

Si 
/(k) 



[ A/2 -70/ (k) \ 
\--Tor (k) -A/2 )' 

( 1 *o/(k)\ 

V *o/* (k) 1 J ' 

,ifc v a/V3 + 2e- lk ^ a/2 ^ cos{k x a/2). 



where <I>p = J e lp r ^(r)d 3 r is the Fourier image of an 
atomic orbital <I> (r) , and the wave vector component par- 



Here, parameter 70 describes the strength of nearest- 
neighbor hopping yielding the Fermi velocity v — 
(a/3/2) aj /h [39 , and a is the lattice constant. The pa- 
rameter so <C 1 describes non-orthogonality of orbitals, 
A = 6a — tB describes a possible asymmetry between A 
and B sites (thus opening a gap |A|). Note that here 
we neglected next-nearest neighbor hops which do not 
produce any visible change in the calculated spectra. 

Figure 1(c) shows constant-energy intensity patterns 
("maps") at 1.45eV below the charge-neutrality point 
|36j in monolayer graphene, plotted as a function of pho- 
toelectron wave vector pii = (j> Xl p y ) parallel to the sur- 
face of graphene, covering the whole Brillouin zone. Since 
the patterns in the vicinity of each Brillouin zone corner 
are the same, but rotated with respect to each other, we 
describe in detail what is happening around one Brillouin 
zone corner. Fig. 2(b) shows a series of plots demonstrat- 
ing the evolution of the constant-energy map with energy 
for the valley K + = (47r/3a,0). Each plot is for a dif- 
ferent fixed energy from the charge-neutrality point with 
energies above (below) on the left (right) hand side. For 
states above the charge-neutrality point (left), the an- 
gular variation is cos 2 (y/2) where <p is the angle of the 
momentum measured from the center of the valley: com- 
parison with Eq. ([3| yields cf> = ip, illustrating that the 
isospin is parallel to the momentum a.iii — 1. Fig. 2(b) 
(right) shows the intensity for emission from states below 
the charge-neutrality point in monolayer graphene. In 
this case, the patterns are flipped with respect to those 
of the left and comparison with Eq. (J3J) yields (j) = ip + tt 
indicating that the isospin is antiparallel to the momen- 
tum cr.ni = — 1. 

Figure 2(c) shows the development of the fixed-energy 
intensity pattern as asymmetry of on-site lattice ener- 
gies A = 6a — £b increases, opening a small gap in the 
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FIG. 2: (a) schematics of the low-energy bands e = ±hvq 
near the .K+ point in the absence of intralayer asymme- 
try A. (b) The intensity of photoemission from states at 
a fixed energy close to the charge-neutrality point in mono- 
layer graphene, plotted as a function of photoelectron wave 
vector qy = (q x ,q y ) parallel to the surface of graphene in the 
vicinity of valley pn = (47r/3a,0). Each plot corresponds to 
a different energy with respect to the charge-neutrality point, 
either above (left hand side) or below (right). Parameter val- 
ues are 70 = 3.0eV, so = 0.129, A = 0, and the energy 
width F varies as the energy divided by six. (c) The develop- 
ment of the intensity pattern for emission at a fixed energy 
O.leV below the charge-neutrality point in the vicinity of val- 
ley P|| = (4-7r/3a, 0) in monolayer graphene as intralayer asym- 
metry increases in magnitude A/2 = 0.05, 0.075, O.leV. Pa- 
rameter values are 70 = 3.0eV, so = 0.129, and the Lorentzian 
energy broadening F — 0.0167. 



smaller and the intensity anisotropy vanishes [Fig. 2(c), 
middle and right]. Thus, the opening of an AB asymme- 
try gap in the monolayer spectrum is accompanied by the 
loss of the chirality-related anisotropy of ARPES angular 
maps at low energies, which can be used not only as a 
strong test for the symmetry-breaking effect induced by 
a SiC substrate [TS] , but also as a probe of wave function 
mixing by AB asymmetry or trigonal warping. Unlike 
AB asymmetry, whose effects on the ARPES intensity 
are strongest near the charge-neutrality point, trigonal 
warping affects the ARPES spectral intensity only at very 
large energies. 

For finite A, an analytical, approximate description 
of the ARPES intensity can be developed as follows. 
For electronic energies much less than the 7r-band width 
(«a< 1), 



/(k) 



V3c 



(q x - iq y ) + -^(q x + nyf , 



(5) 



and 



Til 



A/2 



mW 2 



vn-fi^) 2 -A/2 



7r = hq x + ihq. 



v 



where fi = 7oa 2 /8fi. 2 describes the strength of trigonal 
warping (we assume that \xhq <C v) . This determines the 
spectrum 



e q w s\ h 2 v 2 q 2 — 2^fivh 3 q 3 cos Sip + fi 2 h 4 q 4 



A 2 



(6) 



where s — 1 (s = — 1) stand for the conduction (valence) 
band index, and, using Eq. Q, to the ARPES angular- 
dependent intensity 



1 



hvq 



cos (29) - ^ cos (26 - 3fr>) 



X 5 (E p + A — e q — huj) Sq t 



q,p»-K±-G 



(7) 



where 9 = ^-f (m x - m 2 ) + f (1 - s£). Eqs. (|6 7} con- 



tain the full dependence on valley £ = ±1 and reciprocal 
lattice vector (mi,m 2 ) indices [30]. 



III. PHOTOEMISSION FROM BILAYER 
GRAPHENE 

The use of ARPES to determine the sign of 
interlayer coupling parameter 71 



spectrum. There are two principle effects on the ARPES 
spectrum. For energy gaps \A\/2 -C |e|, the mixing of 
the wave functions between A and B sites destroys the 
perfect cancellation of the ARPES intensity, so that the 
ratio between the maximum and minimum intensity be- 
comes finite [Fig. 2(c), left image]. As the gap increases 
towards the energy probed, the ARPES contour becomes 



Bilayer graphene [THJ H3J consists of two coupled 
hexagonal lattices with inequivalent sites A\, Bl and 
A2, B2 in the first and second graphene sheets, respec- 
tively, arranged according to Bernal (^42-i?l) stacking 
[23] as shown in Fig. 3(a). As in the monolayer, the Bril- 
louin zone has two inequivalent degeneracy points K± 
which determine two valleys centered around zero en- 
ergy in the electron spectrum. Near the center of each 



valley the electron spectrum consists of four branches, 
Fig. 3(b), with two branches describing states on sublat- 
tices A2 and Bl that are split from zero energy by about 
±|7i|, determined by the inter layer coupling 71, whereas 
two low-energy branches are formed by states based upon 
sublattices Al and B2. 

To model bilayer graphene we use a tight-binding 
Hamiltonian matrix H2 and overlap-integral matrix 
§2 that operate in the space of coefficients tp T = 
Wmi,Vb2,Va2,V>Bi) at valley K+ [23l EH E] : 



H 2 = 



S 2 



£A1 

73/ (k) 
74/* (k) 

V -7o/* (k) 
/ 1 



s /(k) 

V s f* (k) 



73/* (k) 

-7o / (k) 
74/ (k) 



1 sof* (k) 
1 

si 



74/ (k) 
■7o/* (k) 

7i 



-7o/(k)\ 
74/* (k) 
7i 



*o/(k)\ 




si 
1 



(8) 



/ 



We adopt the notation of the Slonezewski- Weiss- McClure 
model [2j 13] that is often used to describe bulk graphite, 
in order to parameterize the couplings relevant to bilayer 
graphene |39j . Nearest-neighbor coupling within each 
plane is parameterized by coupling 70 [v — (y3/2\ ajo/h] 
and interlayer A2-BI coupling is described by 71. The 
parameter 73 describes direct AI-B2 interlayer coupling 
which leads to an effective velocity V3 = (v3/2J aj^/h 
representing the magnitude of trigonal warping, partic- 
ularly relevant at low energy [we assume that hvq 3> 
71(1)3/1;)]. Parameter 74 describes A1-A2 and B1-B2 
interlayer hopping. Using 74 = 0.044eV [38j [42] we 
found no noticeable effect of 74 on the ARPES plots, 
and, for simplicity, we use 74 = throughout. Other 
weaker tunneling processes including next-nearest neigh- 
bor hopping are also neglected. The parameter s\ de- 
scribes non-orthogonality terms arising from overlaps be- 
tween orbitals on different layers. Following numerical 
analysis, we found it also has no noticeable effect on the 
ARPES plots In the following angular maps we use 
s% = 0. We note that some works on bilayer graphene 
use different definitions of the tight-binding parameters 
(for example, 73 is defined with an additional minus sign 
in [23J |5HJ |3T] , but it has no effect on their conclusions). 

The Bloch function amplitudes i[> T = 
(i>Ai,ipB2,ipA2,^Bi) and band energy e q , found 
using the Hamiltonian Eq. ([8]), can be used to model 
the photoemission intensity. Figure 3(c) shows constant- 
energy maps at 1.45eV below the charge-neutrality point 
|36j in bilayer graphene, with the plot on the left hand 
side showing values of py covering the whole Brillouin 
zone. The patterns in each valley are the same, but 
rotated with respect to the others, so that we can focus 
on one of them, highlighted in Fig 3(c). The anisotropy 
of the bilayer pattern at this energy is similar to that 
of the monolayer, Fig. 1(c), because the energetic width 
r ~ 1 71 1 obscures features associated with the presence 




FIG. 3: (a) schematic of the bilayer lattice containing four 
sites in the unit cell: ^41 (white circles) and Bl (black) in 
the bottom layer, and A2 (white) and B2 (black) in the top 
layer, (b) Schematic of the hexagonal Brillouin zone with two 
inequivalent valleys K± and the low energy bands in the ab- 
sence of lattice asymmetry. The energy band index a = ±1, 
Eq. (Ill, is shown explicitly for the case 71 < 0. (c) The inten- 
sity of photoemission from states at a constant energy 1.45eV 
below the charge-neutrality point in bilayer graphene, plot- 
ted as a function of photoelectron wave vector py = {p x ,p y ) 
parallel to the surface of graphene for py covering the whole 
Brillouin zone (left) and plotted as a function of photoelectron 
wave vector q = {q x ,q y ) in the vicinity of valley K+ (right) 
[note that the origin and scale of py and q are different]. 
Here we use parameter values 70 = 3.0eV, 71 = — 0.35eV, 
73 = -0.15eV, 74 = O.OeV, s = 0.129, A = U = 0, and the 
energy width T = 0.24eV. 



of two bands. To observe differences between the two 
materials, we need to consider the ARPES patterns at 
energies closer to the charge-neutrality point. 

Figure 4 illustrates the evolution of the intensity pat- 
tern with energy. At energies greater than the interlayer 
coupling, |e q | > I71I (the top two patterns), there are 
two ring-like patterns, each corresponding to photoemis- 
sion from states in two bands, whereas, for low-energies, 
l e g I < l7il (the bottom two patterns), there is a single 
ring corresponding to emission from the degenerate band 
only. Although these plots have been obtained using a 
complete bilayer Hamiltonian Eq. (JsJ> , it is convenient 
to discuss salient features of the results in Fig. 4 using 
an analytic formula, obtained by performing a linear- 
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FIG. 4: The intensity of photoemission from states at a 
fixed energy close to the charge-neutrality point in bilayer 
graphene, plotted as a function of photoelectron wave vector 
q = (q x ,q y ) parallel to the surface of graphene in the vicin- 
ity of valley pn = (47r/3a,0). Each plot corresponds to a 
different energy with respect to the charge-neutrality point, 
either above (left) or below (right). Parameter values are 
70 = 3.0eV, 71 = -0.35eV, 73 = -0.15eV, 74 = O.OeV, 
so = 0.129, A — U = 0, and the energy width V varies as 
the energy divided by six. 



in- momentum expansion of /(k), Eq. ([5]), and neglect- 
ing trigonal warping due to AI-B2 interlayer coupling 
(73 = 0), A1-A2 and B1-B2 interlayer coupling (74 = 0), 
and non-orthogonality of orbitals (s± = sq = 0). In this 
case, the four bands in the bilayer spectrum are described 

by 



1 7i I 



1 + 4/i 2 vV/7i 2 + b 



(9) 



where the parameters 
b = ±l 



s = ±1 



identify the four bands: 6=1 for the split bands with 
energy \e q \ > and b = — 1 for the low-energy 'de- 
generate' bands that touch at zero energy, while s = 1 
(s = —1) indicates the conduction (valence) bands. 
Then, the contribution of a given band is 



I -[i + { e q /nv q f] 5 ^ + A 



where 

9(<P) z 



^w) £q,p,|-K ± -Gj 



(10) 



1 + acos(2ip) + 5 at x 



2 



and 



S&71/I71I. 



(11) 



As the value of a, Eq. (Ill, depends on the sign of the 



tight- binding parameter 71, comparison of the angular 
dependence of g{ip) with experimental data provides a 
method to determine the sign of 71 [33]. To demon- 
strate this, we make a comparison with our numerical 
data, plotted in Fig. 4. In this illustration, we assume 
that 71 < 0, which is a natural choice given the z — > — z 
asymmetry of the p z orbitals of carbon. It shows how, 
for this choice of the sign of 71 , the anisotropy of pho- 
toemission angular-maps differ in the split bands and de- 
generate bands at energies above e > and below e < 
the charge neutrality point. Note that changing the sign 
of 71 to positive would lead to an interchange of plots 
illustrating the ARPES behavior at e > and e < 0. 

The most pronounced feature of the ARPES angular 
maps, depicted for 71 < in Fig. 4, is that, for en- 
ergies e > (left side of Fig. 4), photoemission spec- 
tra are dominated by states in the degenerate bands, 
b = — 1, which are nicely described by the intensity profile 
I cx cos 2 ip/2. In contrast, for e < (valence bands, right 
side of Fig. 4), ARPES intensity from the degenerate 
band b = — 1 is weak, whereas the split band, at energies 
e < — 1 7! |, produces a bright, dominant signal. If ob- 
served experimentally, such a behavior of ARPES maps 
in the conduction and valence bands would be indicative 
of a negative sign of the interlayer coupling 71 [39]. If 
the experimentally-observed constant-energy maps were 
interchanged for negative and positive energies, it would 
be evidence for 71 > 0. Although the sign of 71 has di- 
rectly observable consequences for the ARPES pattern, 
tight-binding parameters for graphite published so far 
have assumed 71 > ([43 and references therein). 



B. Electron chirality in the ARPES of bilayer 

graphene and the use of trigonal warping to 
determine the interlayer coupling parameter 73 

The behavior of low-energy particles in bilayer 
graphene is perhaps even more remarkable [10] [13] [23] 
than in a monolayer. The low-energy bands (at energy 
|e| <C |7i|) have a parabolic energy- versus-momentum re- 
lation and they support eigenstates of an operator a.n 2 
with cr.n.2 = 1 for electrons in the conduction band and 
(T.n 2 = —1 for electrons in the valence band, where 
n 2(q) = (cos(2y>), sin(2y>)), which means that they are 
chiral, but with a degree of chirality different from that 
in the monolayer, with the isospin linked to, but turn- 
ing twice as quickly as, the direction of momentum. An 
interpretation of the ARPES constant-energy maps in 
terms of two-source interference Eq. ^ predicts an an- 
gular variation like cos 2 (y) for states above the charge- 
neutrality point and cos 2 (c/? + tt/2) for states below (for 
7i < 0). 

Figure 5 shows the calculated intensity of the photoe- 
mission in the vicinity of valley = (Air /3a, 0) from 
states very close to the charge-neutrality point in bi- 
layer graphene, at energy 0.03eV above (left-hand side) 
and 0.03eV below (right). We consider two different 
signs of the A2-B1 interlayer coupling strength, ji, with 
71 = -0.35eV (top) and 71 = +0.35eV (bottom). For 
71 < and emission from the valence band (top right), 
the pattern is like cos 2 (<^ + tt/2) as expected for two- 
source interference of chiral electron in bilayer graphene. 
As shown in Eq. (10) [and explained in detail in Sec- 
tion |IV| , the intensity from this band is not affected by 
corrections due to the presence of dimcr A2 — Bl or- 
bitals (it has a = —1). For emission from the conduction 
band (top left side of Fig. 5), the interference pattern has 
two peaks, but one of the peaks has about three times 
stronger maximum intensity than the other because of 
the presence of the contribution from dimer A2 — Bl or- 
bitals (this band has a — +1) [H]. The bottom-left and 
bottom-right plots in Fig. 5 show the constant-energy 
maps for 7! > for emission above and below the charge- 
neutrality point, respectively. In this case, the intensity 
pattern for emission from the conduction band (bottom 
left) has two peaks with the same maximum intensity, 
arising from the interference of waves from the Al and 
B2 sublattices. For emission from the valence band (bot- 
tom right) the peaks have different maximum intensities, 
owing to the interference of waves from four sublattices. 

We note that, once the sign of is known, the sign of 
AX-B2 interlayer coupling 73 may also be deduced from 
the orientation of trigonal warping of the intensity pat- 
terns near the charge- neutrality point |39j . In bilayer 
graphene, there are two principle causes of trigonal warp- 
ing. The first is the presence of AX-B2 interlayer coupling 
73 that will tend to dominate at low energy, the sec- 
ond is higher-in-momentum terms in the function /(k) 
that will be important at large energy. The latter causes 
trigonal warping in monolayer graphene while the former 
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FIG. 5: The intensity of photo-emission in bilayer graphene 
for fixed energy very close to the charge-neutrality point in 
the vicinity of valley pn = (47r/3a,0): for states with en- 
ergy 0.03eV above the charge-neutrality point (left) and states 
with energy 0.03eV below the charge-neutrality point (right). 
We consider different signs of the A2-B1 interlayer coupling 
strength, 71, with 71 = — 0.35eV (top) and 71 = +0.35eV 
(bottom). Other parameter values are 70 = 3.0eV, 73 = 
-0.15eV, 74 = O.OeV, F = 0.005eV, and s = 0.129 [37]. 



is not present in a monolayer. At large energies, when 
the higher-in-momentum terms dominate, the orientation 
of trigonal warping is the same in bilayer graphene [e.g. 
Fig. 3(c)] as in a monolayer [e.g. Fig. 1(c)] whereas, at 
low energy, the orientation of trigonal warping in a bi- 
layer depends on the sign of parameter 73 (assuming that 
the sign of 71 is known) . 

The orientation of trigonal warping flips on changing 
the sign of 71 as seen by comparing the top and bottom 
plots in Fig. 5. At very low energy, e q ,hvq <C 71, and in 
the absence of lattice asymmetry, the energy eigenvalues 
|2"B1 are 



±,KV-2C ^ W co S 3^+^y. 
V 7i 7i 



(12) 



where W3 = — (v3/2) ajs/h. This expression illustrates 
that the angular dependent factor, producing trigonal 
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A = 0, U = 



A*0,U = 



A = 0, U * 
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FIG. 6: The band structure of bilayer graphene in the vicinity 
of a valley for no lattice asymmetry (left), substrate-induced 
asymmetry A = 71 < (U = 0) (center), interlayer asymme- 
try U = 7i/2 (A = 0) (right). For clarity, we use large values 
of asymmetry. The energy band index a = ±1, Eq. (Ill, is 



shown explicitly for the case 71 < on the left hand side. 



warping, depends on the sign of the ratio 73/71. In this 
paper we usually choose ji < and 73 < to illustrate 
the possibility that the orientation of trigonal warping is 
different at lower energies [e.g. Fig. 5 (top)] than that at 
higher energies. 



Substrate-induced asymmetry in bilayer 
graphene 



The Hamiltonian H2, Eq. (|8|, takes into account the 
possibility of different on-site energies through its diag- 
onal components. Their effect may be understood by 
considering the eigenenergies exactly at the center of the 
valley where /(k) = 0, namely e = €ai, e = £52, or 



e = I (eA2 + (bi) ± V I ( e A2 - em 



Below, we distinguish between two types of asymmetry 
in bilayer graphene parameterized using A = cai — £bi 
the difference between on-site energies of adjacent atoms 
in the bottom layer due to the presence of a substrate, 
and interlayer asymmetry U — [(cai + cbi) — (eA2 + 
£B2)]/2 between on-site energies in the two layers arising 
from a doping effect and charge transfer to the substrate 

In Fig. [6] (center), the band structure in the vicinity 
of the K point is plotted in the presence of substrate- 
induced asymmetry A = e^i — esi [the plot is shown for 
71 < and A < 0] . This type of asymmetry introduces 
a gap ~ I A |/2 as well as electron-hole asymmetry. In 
Fig. [6] (right), the band structure in the vicinity of the K 
point is plotted in the presence of interlayer asymmetry 
U = [(e A i + est) - (e^2 + es 2 )]/2. It does not break 
electron-hole symmetry, but introduces a gap ~ |Z7|. 

As illustrated in Fig. 7, the constant energy maps are 
sensitive both to the magnitude and sign of the asymme- 
try A. The plots on the left (right) side of Fig. 7 show 
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FIG. 7: Left (right) shows the development of the intensity 
pattern for emission at a fixed energy O.leV above (below) 
the midgap energy in bilayer graphene in the vicinity of val- 
ley P|| = (4n/3a, 0) as intralayer asymmetry | A| increases [the 
energy with respect to the charge-neutrality point is also indi- 
cated]. The plots show A = (top), negative A/2 = — 0.15eV 
(middle), and positive A/2 = 0.15eV (bottom). Parame- 
ter values are 70 = 3.0eV, 71 = — 0.35eV, 73 = — 0.15eV, 
74 = O.OeV, s = 0.129, V = 0.0167eV. 



constant energy maps for photoemission from conduction 
(valence) band states at energy O.leV above (below) the 
midgap energy. The top two plots are for no asymmetry 
A = 0, the middle two plots show negative asymme- 
try A/2 = — 0.15cV, and the bottom two show positive 
asymmetry A/2 = +0.15eV. As for the monolayer, one 
effect of asymmetry A is to impair the two-source inter- 
ference resulting in a weakening of the angular anisotropy 
of the intensity pattern. The 'Mexican hat' structure of 
the valence band for negative A, and the conduction band 
for positive A, is manifested in the larger ARPES con- 
tour for emission from these states (shown on the mid- 
dle right and the bottom left, respectively) as opposed 
to their counterparts in the other band (middle left and 
the bottom right, respectively). Experimentally, such a 
difference in the size and nature of the ARPES contour 
for emission from conduction or valence bands (at the 
same distance from the midgap energy) would indicate 
the presence and sign of intralayer asymmetry A. 

Since interlayer asymmetry U results in a gap ~ \U\ 
that preserves electron-hole symmetry and does not de- 
pend on the sign of U, the constant-energy photoemission 
maps are sensitive to the magnitude of U but not its sign. 
Figure 8 shows constant energy maps for photoemission 
from conduction band states at energy O.leV above the 
charge-neutrality point (left hand side) and from valence 
band states at energy O.leV below the charge-neutrality 
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FIG. 8: The development of the intensity pattern in bilayer 
graphene in the vicinity of valley py = (4-7r/3a,0) as inter- 
layer asymmetry (7/2 = 0, 0.075, 0.1, 0.125, 0.15eV increases, 
for emission from states with energy O.leV above the charge- 
neutrality point on the left hand side, and energy O.leV below 
the charge-neutrality point on the right. Parameter values 
are 70 = 3.0eV, 71 = -0.35eV, 73 = -0.15eV, 74 = O.OeV, 
s = 0.129, T = 0.0167eV. 



point (right) as U increases in magnitude. Generally, 
the effect of asymmetry U is to weaken the angular 
anisotropy of the intensity pattern. Both the conduction 
and valence band in the vicinity of the charge neutrality 
point display a 'Mexican hat' structure 23, 24, 25 1, lead- 
ing to relatively large ARPES contours (and additional 
features as seen in the plot at U/2 = O.leV), in con- 



trast to asymmetry A, Fig. 7, where the 'Mexican hat' 
structure appears in one band only. 



IV. DISCUSSION 

Our numerical data for bilayer graphene, Figs. 4,5, 
show that the anisotropy of photoemission angular-maps 
differ in the split bands and degenerate bands at energies 
above e > (left side of Fig. 4) and below e < (right 
side of Fig. 4) the charge neutrality point. These plots 
may be interpreted in terms of the interference of photo- 
electron waves emitted from four non-equivalent sublat- 
tices. In fact, for two of the four bands, the parameter 
a = —1 in Eq. 
on the 'dimcr' 



(10 1, so that the contribution of orbitals 
sites A2 and Bl cancel, leaving only the 



contribution of two terms e v in g((p) that arise from 
orbitals on sites Al and B2. For the other two bands, 
a = 1, the contribution of orbitals on the 'dimer' sites 
A2 and Bl to g(ip) do not cancel, but interfere with the 
contribution of orbitals on sites Al and B2, producing a 
different angular dependence and greater peak intensity 
than for a = — 1. 

As the value of a, Eq. (Ill, depends on the sign of 71, 



comparison of the angular dependence of g(f) with ex- 
perimental data provides a method to determine the sign 
of 71 [39 . To demonstrate this, we make a comparison 
with our numerical data, plotted in Fig. 4. For the sign 
of 71 that we adopt in the numerics (71 < 0), the split 
band above the charge- neutrality point has a = — 1 so 
the intensity from this band appears as a very faint ring 
(that of smaller radius) in the plot at energy e = +0.5eV 
in Fig. 4. The degenerate band at this energy, however, 
has a = 1 so the intensity from it appears as the ring of 
larger radius with larger peak intensity. As the energy 
drops below |7i|, left side of Fig. 4, the contribution of 
the split band disappears to leave only the ring arising 
from the degenerate band with a = 1. The energy has to 
approach the charge-neutrality point before the contribu- 
tion of the dimer sites A2 and Bl, small in the parameter 
e q /hvq ~ v/eg/71, weakens to reveal an anisotropy pat- 
tern characteristic of two source interference in bilayer 



graphene, as explained in Section III B [44 . 

The picture is quite different for energies below the 
charge-neutrality point (right side of Fig. 4) . In this case 
the split band has a = 1 so that the intensity from it 
appears as the ring (of smaller radius) with larger peak 
intensity at energy e = — 0.5eV in Fig. 4. The degenerate 
band has a = — 1 so the intensity from it appears as the 
fainter ring (that of larger radius) at energy e = — 0.5eV. 
As the energy increases above — I71I, right side of Fig. 4, 
the contribution of the split band disappears to leave only 
the ring arising from degenerate band with a = — 1. This 
is why the intensity pattern ~ cos 2 (<y9) is much easier to 
detect below the charge-neutrality point than above it. 
In fact, whether it is easily visible above or below the 
charge-neutrality point depends on the sign of 71 (here 
we chose 71 < 0) so the experimental observation of the 
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FIG. 9: The development of the intensity pattern in bilayer 
graphene for emission from states with energy O.leV below 
the charge-neutrality point in the vicinity of valley pn = 
(47r/3a, 0) as attenuation [described by factor exp(— 2z + 2i/3) 
where (z,/3) are real parameters] of waves from the bot- 
tom layer, increases. Parameter values are 70 = 3.0eV, 
71 = -0.35eV, 73 = -0.15eV, 74 = O.OeV, s = 0.129, 
T = 0.0167eV. 

anisotropy ~ cos 2 (ip) will provide a way to determine the 
sign of A2-BI intcrlayer coupling 7! in bilayer graphene. 

Finally, we note that the anisotropy of the constant- 
energy maps may be influenced by other factors not mod- 
eled here. For example, when the component of photo- 
electron momentum perpendicular to the bilayer sample 
p z is large, we expect that photoelectron waves emitted 
from the bottom layer will be attenuated with respect to 
those emitted from the top layer. To obtain an impres- 
sion of the typical kind of effect, we introduced an expo- 
nential attenuation [described by factor exp(— 2z + 2i/3) 
where (z, [3) are real parameters] of waves from the bot- 
tom layer. As shown in Fig. 9 for photoemission from 
states at energy O.leV below the charge-neutrality point 
in bilayer graphene, the attenuation results in a destruc- 
tion of the double-peaked intensity pattern and the phase 
factor (3 has the effect of rotating the whole pattern. 

V. CONCLUSIONS 

Using Fermi's Golden Rule, we modeled the anisotropy 
of the intensity of photoemission constant-energy maps 
at low energy in graphene and demonstrated that the 
anisotropy is a manifestation of electronic chirality. In 
monolayer graphene, photoemission may be viewed as a 
two-source interference experiment, a la Young's double 
slits, the sources being two inequivalcnt lattice sites in 
the unit cell. The resulting intensity ~ cos 2 (</?/2) dis- 
plays a single-peaked dependence on the direction of mo- 
mentum described by angle <p. In bilayer graphene, the 
interference of emitted photoelectron waves from four 
atomic sites produces single- or double-peaked constant- 
energy maps, depending on the energy of the initial 



state in graphene. The marked contrast between the 
anisotropy for emission from the conduction or the va- 
lence band at energies below the A2-B1 intcrlayer cou- 
pling strength, parameterized by 71 [39 , provides an ex- 
perimental method to determine the magnitude and sign 
of parameter 71. 

The shape of the photoemission constant-energy maps 
is determined by the trigonal warping effect in graphene. 
In monolayers and bilayers, the iso-energetic line changes 
from an almost-circular to a triangularly-warped shape 
as the energy increases: the extent of such warping is 
controlled by the dimcnsionlcss parameter qa, where q is 
the magnitude of the wave vector measured from the cen- 
ter of the valley and a is the lattice constant. In bilayer 
graphene, strong trigonal warping may also occur at low 
energy because of A\-B2 interlayer coupling, parameter- 
ized by 73 [32] , and the observation of this latter trigonal 
warping provides an experimental method to determine 
the magnitude and sign of parameter 73. 

Measurement of the anisotropy of the intensity of pho- 
toemission constant-energy maps provides a method to 
characterize realistic graphene samples. As an exam- 
ple, we take into account substrate-induced asymme- 
try that impairs the two-source interference in mono- 
layer graphene, resulting in a weakening of the angular 
anisotropy of the intensity pattern. Analysis of recent 
experimental data |15l I16j in terms of the anisotropy of 
constant-energy maps may help to shed light on the pos- 
sible presence of asymmetry in graphene grown epitaxi- 
ally on SiC substrate. In bilayers, both substrate-induced 
asymmetry and interlayer asymmetry alter the interfer- 
ence pattern: we describe measurable differences between 
them. This illustrates the potential of photoemission in 
the future characterization of few-layer graphene sam- 
ples. 
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